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2.5 Stokes flow past a sphere 

[Refs] 

Lamb: Hydrodynamics 

Acheson : Elementary Fluid Dynamics, p. 223 ff 

One of the fundamental results in low Reynolds hydrodynamics is the Stokes solution for 
steady flow past a small sphere. The apUicatiuon range widely form the determination of 
electron charges to the physics of aerosols. 

The continuity equation reads 

V-q = (2.5.1) 
With inertia neglected, the approximate momentum equation is 

= -— + I/VV (2.5.2) 

P 

Physically, the presssure gradient drives the flow by overcoming viscous resistence, but does 
affect the fluid inertia signiflcantly. 

Refering to Figure 2.5 for the spherical coordinate system [r,6,(f)). Let the ambient 
velocity be upward and along the polar {z) axis: {u,v,w) — (0, 0, H^). Axial symmetry 
demands 

d 

— = 0, and q= {qr{r,9),q0{r,9),Q) 

Eq. (2.5.1) becomes 

As in the case of rectangular coordinates, we define the stream function tjj to satisify the 
continuity equation (2.5.3) identically 

1 dij 1 d^ip 

= 2 • /j ^' ^0 = (2.5.4) 

sm ou r sm or 

At infinity, the uniform velocity W along z axis can be decomposed into radial and polar 
components 

1 dip 1 dib , , 

q, = Wcose = qe = -W sine ^ ^ ~ 2.5.5 

sm oB r sm or 



Figure 2.5.1: The spherical coordinates 



The corresponding stream function at infinity follows by integration 

W 



ih — —r^ sin^ 9, r ~ oo 
^ 2 



Using the vector identity 



and (2.5.1), we get 



V X (V X = V(V-g) - VV 



VV= -V X (V X = -V X C 



Taking the curl of (2.5.2) and using (2.5.8) we get 

V X (V X C) = 

After some straightforward algebra given in the Appendix, we can show that 



(2.5.6) 

(2.5.7) 
(2.5.8) 

(2.5.9) 



g = V X 



r sin 9 



(2.5.10) 



and 



(=Vxg=VxVx 



^60 



d'^ip ^sm.9 d ( 1 dil)^ 



rsin9j r sin ^ \ dr"^ 39 \ sin 9 89 



(2.5.11) 



Now from (2.5.9) 

VxVx(Vxg) = VxVx 



V X V X 



rsin^ 
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hence, the momentum equation (2.5.9) becomes a scalar equation for ■0. 



' sin^ (9/1 



V' = 



The boundary conditions on the sphere are 

= qe = ^ on r = a 

The boundary conditions at oo is 

W 



sin^ e 



(2.5.12) 
(2.5.13) 
(2.5.14) 
(2.5.15) 



Let us try a solution of the form: 

^{r,e) = /(r)sin2^ 
then / is governed by the equi-dimensional differential equation: 

whose solutions are of the form /(r) oc r", It is easy to verify that n=— 1,1, 2, 4 so that 

A 



/ = 



(2.5.16) 



/(r) = - + Br + Cr^ + Dr^ 
r 



or 



■0 = sin^ 6 



- + Br + Cr^ + Dr^ 
r 



To satisfy (2.5.14) we set D = 0, C = W/2. To satisfy (2.5.13) we use (2.5.4) to get 

W_ 

T 



W A B A B 

gr = = — + ^ + - = 0, qg = Q = W- — + - = 



a"^ a 



a"^ a 



Hence 



A = -Wa\ 
4 



B = --Wa 
4 



Finally the stream function is 



W 



3ar 



r" H 

2r 2 



sin^^ 



(2.5.17) 



Inside the parentheses, the first term corresponds to the uniform flow, and the second term 
to the doublet; together they represent an inviscid flow past a sphere. The third term is 
called the Stokeslet, representing the viscous correction. 
The velocity components in the fluid are: (cf. (2.5.4) : 

3a 



= Videos 6' 1 + 
qQ — — VFsin^ 1 — 



2r3 2r 
3a 

4^3 4^ 



(2.5.18) 
(2.5.19) 
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2.5.1 Physical Deductions 

1. Streamlines: With respect to the the equator along 9 = 7i/2, cos 6' and Qr are odd while 
sin 9 and qo are even. Hence the streamlines (velocity vectors) are symmetric fore and 
aft. 

2. Vorticity: 

(ld{rqg) 3 sin6'^ 



3. Pressure : Prom the r-component of momentum equation 

^ /f^^^g^^^ y ^ ^y ^ 

or r"* 

Integrating with respect to r from r to oo, we get 

3/^a 
2~ 



P = Poo- T^^^cos^ (2.5.20) 



4. Stresses and strains: 

On the sphere, r — a, Srr — hence cr^r = and 



On the other hand 



Hence at r = a: 



— p + = — Poo + cos 9 (2.5.21) 

2 Qj 



= r—('l2\+-^ = --^sm9 
dr \ r J r 89 2 r'^ 



Tre = cTre = l^'^re = -^^^sin6l (2.5.22) 



The resultant stress on the sphere is parallel to the z axis. 

Lj; = Trr cos 9 — TrQ Sm 9 = — Poo COS + 

2 a 

The constant part exerts a net drag in z direction 

r27r rn 3 nW 



D= ad(p d9sm9E,^^ -^Atto^ = 67riJ,Wa (2.5.23) 

Jo Jo 2 CL 



This is the celebrated Stokes formula. 
A drag coefficient can be defined as 



r, D 67r/.M^» 24 

~ ^pWW ~ \pWW ~ ~ Red ^ ' 
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5. Fall velocity of a particle through a fluid. Equating the drag and the buoyant weight 
of the eparticle 

4:71 

GTTiiWoa = —a^{Ps - Pf)g 



hence 



^ \^ Pf ) \^ Pf J 



in cgs units. For a sand grain in water, 

Ap _ 2.5 - 1 _ 
Pf 1 



1.5, 1/ = 10-^cm7s 



Wo = 32, 670 a^cm/s 
To have some quantitative ideas, let us consider two sand of two sizes : 

a = lO^^cm = lO^^m : Wo = 3.27cm/s; 

a = IQ-^cm = 10"^ = 10//m, Wo = 0.0327cm/s = 117cm/hr 



(2.5.25) 



For a water droplet in air, 

Ap 1 



Pf 10-^ 



10^ = 0.15 cmVsec 



then 



Wo^ 



(217.8)10^ 2 
015 



in cgs units. If a = 10 ^ cm = 10pm, then M/q = 1.452 cm/sec. 



Details of derivation 



Details of (2.5.10). 



(2.5.26) 



V X 



1 



r sin 9 I sin 9 



Cr eg r sm fc 

d_ d_ A. 

dr ae d<t> 

V 



' 1 9V' ^ 

sin 9 d9 



' 1 di>' 
r sin ^ ()r 
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Details of (2.5.11). 



V X V X 



r sin 9 



V xq 



1 



sin 9 



Cr reo 
A. A. 

dr 89 
1 §2^ -1 dji 
sin 6 dr 



rsin9 



d'^(f) sin 6* d 



Qj.2 



+ 



r sin 9e(^ 

d_ 

d<t> 





r2 d9 \ sin 9 89 



